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We investigate phase transitions in scalar field theories using the functional renormalization group 
(RG) equation. We analyze a system with U(2) x U(2) symmetry, in which there is a parameter 
A2 that controls the strength of the first-order phase transition driven by fluctuations. In the limit 
of A2 — > 0, the U(2) x U(2) theory is reduced to an 0(8) scalar theory that exhibits a second- 
order phase transition in three dimensions. We develop a new insight for the understanding of the 
fluctuation-induced first-order phase transition as a smooth continuation from the standard RG flow 
in the 0(8) system. In our view from the RG flow diagram on coupling parameter space, the region 
that favors the first-order transition emerges from the unphysical region to the physical one as A2 
increases from zero. We give this interpretation based on the Taylor expansion of the functional RG 
equations up to the fourth order in terms of the field, which encompasses the e-expansion results. 
We compare results from the expansion and from the full numerical calculation and find that the 
fourth-order expansion is only of qualitative use and that the sixth-order expansion improves the 
quantitative agreement. 
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I. INTRODUCTION 

Phase transitions occur in many different physical sys- 
tems. In this paper we shall address a fluctuation- 
induced first-order phase transition or the Coleman- 
Weinberg potential [ll-fllj in the system with U(2) x U(2) 
symmetry, as well as a second-order one in the system 
with 0(8) symmetry (l2l - [24| as a particular limit of the 
U(2) x U(2) theory. This particular choice of the global 
symmetries is motivated by those relevant for phase tran- 
sitions in the strong interactions. 

The fundamental theory of the strong interactions is 
Quantum Chromodynamics (QCD), which describes the 
dynamics of quarks and gluons. It is known that at a tem- 
perature T of the same order as the typical QCD scale, 
Aqcdi QCD undergoes two types of phase transitions — 
the chiral pha se transition and the quark deconfinement 
transition [25| . In particular at the chiral phase transi- 
tion the behavior of the system should be characterized 
by global chiral symmetry @, [26| . 

For iVf flavors of massless quarks the QCD Lagrangian 
is invariant under a flavor rotation of the U(iVf)L x 
U(iVf) R symmetry for left- and right-handed quarks. Be- 
cause of the axial anomaly, the subgroup U(1)a is ex- 
plicitly broken (down to Z2Ar f ) on the quantum level, so 
that in this case the (continuous) symmetry is reduced to 
SU(A^f) L x SU(AT f ) R x U(l) v - Here, the vector symmetry 
U(l)v corresponds to the conserved baryon charge and 
is not broken in the normal phase. 

In the QCD vacuum chiral symmetry is spontaneously 
broken by a nonzero chiral condensate (27j according to 



SU(AT f ) L x SU(iV f ) R -> SU(AT f ) v . In the two-flavor case 
the symmetry breaking pattern, SU(2)l x SU(2) r — > 
SU(2)y, is equivalent to SO(4) — > SO(3), and hence can 
be mapped onto a (scalar) meson model with 0(4) sym- 
metry [28f . This provides a motivation to investigate 
phase transitions of the O(N) scalar theory in the finitc- 
T formalism in d — 4 dimensions, or in the zero-T formal- 
ism in d — 3 dimensions, if we assume that the critical 
behavior is described with a dimensionally reduced the- 
ory [2^]. In this paper we address the Renormalization 
Group (RG) flow of the theory in the dimensionally re- 
duced description only. This can be justified as follows: 
At finite temperature, the Euclidean time direction in 
a d-dimcnsional path-integral description becomes com- 
pactificd. Because of this finite extent in the Euclidean 
time direction, in a region close enough to a critical point, 
the dominant long-range fluctuations will be unable to 
propagate in the time direction and the system will in ef- 
fect become (d— 1) dimensional. For a finite-T theory in 
d = 4 space-time it is therefore expected that the critical 
behavior can be described in terms of an effective d = 3 
dimensional theory at T = with its couplings (which 
are T-dependent) set close to their critical values. 

In contrast to the situation described above, if the axial 
symmetry is effectively restored close to the transition 
temperature (e.g . because of the instanton suppression 
at high T [30l - [32| ). the relevant symmetry for the Nf = 2 
transition is then U(2)l x U(2) r , which leads us to the 
analysis on the U(2) x U(2) scalar theory. 

Assuming that a particular phase transition is of sec- 
ond order and displays critical behavior, one can rely on 
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universality arguments to describe the system in terms of 
an effective Landau-Ginzburg- Wilson functional which is 
determined solely from the symmetries and the dimen- 
sionality of the system. The description in terms of such 
an effective theory is often much simpler than the origi- 
nal microscopic theories, but still captures the long-range 
phenomena at the critical point. 

In the investigation of critical phenomena, the correct 
treatment of low-momentum modes is of paramount im- 
portance. Close to a critical point, physical properties of 
a system are dominated by these low-lying soft modes. 
Since the system has no intrinsic length scale at the crit- 
ical point, phenomena on all length scales contribute and 
non-trivial anomalous dimensions can appear. In the case 
of scalar theories, the anomalous dimension at the criti- 
cal point is generally small, which can justify a simplified 
treatment. 

Rcnormalization group methods have been used to deal 
with the infrared (IR) divergences which appear in the 
vicinity of critical points [H[. Historically speaking, the 
e-expansion around d = 4 dimensions for theories with 
upper critical dimension d = 4 has been of great im- 
portance [13, H|[. While the RG methods are able to 
deal with IR divergences, calculations of this type are 
still pcrturbative and rely on an additional expansion in 
terms of small coupling constants, which can be justified 
by the fact that the coupling constants are of order of 
e in the vicinity of the fixed points. There is a count- 
less number of works on the application of RG methods 
to the second-order phase transition in the O(N) theory 

On the other hand, the theory with U(2) x U(2)- 
symmetry requires more investigation. In the context 
of the chiral phase transition in QCD, such an effective 
theory has been extensively discussed by Pisarski and 
Wilczek 0, @, [H[ ■ In their results the e-expansion anal- 
ysis does not find an IR stable fixed point which could 
control a continuous transition. This indicates that the 
transition is presumably driven by fluctuations to be of 
first order. [See also [!, H, Hfil]]- 

This question of the order of the chiral phase transition 
including the possibility of the effective U(1)a symmetry 
restoration has been further analyzed in and beyond the 
e-expansion in e.g. [37T - |39| using field-theoretical methods 



and in |40l - |43j using functional RG methods 
well as the finite-T lattice-QCD simulation [47]. 
related first-order phase-transitions have been investi- 
gated in [9MTTI HU with the functional RG. In a gauge 
theory setting, the functional RG has been used to esti- 
mate effects of the anomaly on chiral symmetry break- 
ing [4!| and on the phase transition line [5(| in QCD. 
Also, anomal y e ffects on the mass spectrum have been 
described in [5l| by means of non-perturbative Dyson- 
Schwinger equations. 

The functional RG approach does not rely on an ex- 
pansion in terms of a small coupling and is capable of 
describing the full effective potential. This enables us to 
directly analyze the critical behavior associated with not 



only the second-order but also first-orderphase transi- 
tions. For reviews of the method see e.g. [45l. l46l [52| - [60| . 

We use a functional RG equation to describe the 
behavior of the U(2) x U(2) matrix model at the 
fluctuation-induced first-order phase transition in the 
same way as done by Berges and Wetterich [42j]. We 
establish the connection between the functional RG flow 
equations in the formulation according to [44| (i.e. the 
Wetterich equation) and the e-expansion. While there is 
some literature on the connection between the functional 
RG and pcrturbative calculations 

[mum til, tws con- 
nection we address is a useful starting point to explore the 
general RG flow diagram. Moreover, to our knowledge, 
nobody has explicitly retrieved the e-expansion results 
from the Wetterich equation for the U(2) x U(2) system. 
[See also [H, [6^] . The connection of the e-cxpansion to 
functional RG equations in a related formalism and the 
emergence of the Wilson-Fisher fixed-point in this ap- 
proach have been pointed out in [60}.] 

Since we are motivated by an interest in the QCD 
phase transition at finite temperature, we will adhere 
to the QCD terminology throughout this paper. That 
is, we will refer to the 8-component fields in the 0(8) 
and U(2) x U(2) theory as the cr-ficld for the condensing 
mode which acquires a non-zero expectation value and as 
the 7?-fields for the 3 pseudo-scalar components that are 
Nambu-Goldstone bosons. These 4 mesons have parity 
partners, namely, 77-field and ao-fields. Apart from this, 
nevertheless, our treatment is generic and is completely 
parallel to the earlier works by Berges and Wetterich [421 
and by Berges, Tetradis, and Wetterich [4!|. In the next 
section, as a convenience to readers, we will briefly sum- 
marize our aims in this paper and discuss where they 
advance beyond those earlier works. 

We intend the current calculation as a starting point 
for further investigations of the transition in the frame- 
work of more realistic phenomcnological models for QCD 
thermodynamics [6#T72| . In such models gauge degrees 
of freedom can also be taken into account, but the pos- 
sibility of a fluctuation-induced change of the order of 
the transition has not yet been fully investigated due to 
technical complexity. The inclusion of such effects be- 
yond the mean-field level in these models would be a 
very interesting problem in the future. 

We will organize this paper as follows: We first present 
a summary of our central results in Sec. |TTJ Then we 
make a quick review over the functional RG formalism 
in Sec. Mil We discuss the analytical and numerical cal- 
culations in great details in Sec. IIV1 Section IVl contains 
our concluding remarks. 



II. CENTRAL RESULTS 

The phase transitions of the U(2) x U(2) scalar theory 
have been extensively studied in earlier works in the func- 
tional RG formalism [42|, EH . Although there are minor 
differences in the technical setup, we have performed nu- 
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merical calculations at the same level of approximation 
and truncation as used by Berges, Tetradis, and Wet- 
terich in [H, [43| ■ Before we discuss the details of the 
calculation, we therefore wish to point out the qualita- 
tively new aspects of this work. 

First of all, we would like to emphasize that we use a 
very simple form of the flow equations (given by Eqs. (|39[) 
and (|43|) below) in the local potential approximation 
which does no longer involve momentum integrations. 
This difference from the equations in [42|, |43| arises from 
a different choice of the IR regulator function. Apart 
from the wave-function renormalization, which we ne- 
glect, the information contained in the RG flow equation 
is equivalent. Since the content is equivalent, the sim- 
pler expressions are advantageous, both for an analysis 
of the structure of the RG flow and for numerical calcu- 
lations. In fact, as we will demonstrate explicitly later 
in the discussion, we can reproduce the result from the 
e-expansion in a short calculation of only a few lines and 
make the connection between the functional RG equation 
and the perturbative method very transparent. This is 
certainly not a qualitative breakthrough, but such a tech- 
nical improvement provides useful progress for practical 
calculations. 

Second, in this work we propose a new point of view for 
understanding the fluctuation-induced first-order phase 
transition. If the phase transition is of second order, the 
flow diagram and the fixed-point structure can be studied 
very well in the space of the coupling parameters. For 
instance, in the case of the O(N) scalar theory in three 
dimensions with a field <fi = (<f>i, . . . 4>n) T , the relevant 
operators in terms of the squared field ip — 4> 2 are ip and 
tp 2 . The (IR-scale dependent) effective potential can thus 
be written as 

£4(</?) = ^a4</? + ^W 2 . (l) 

The relevant coupling parameters, /i 2 and Aifc, run as 
the IR scale k becomes smaller and fluctuations are in- 
tegrated out. Using dimensionless variables p, 2 = k~ 2 /j, 2 
and Aife = fc^Aifc, we obtain Fig. [1] Starting from a 
set of differential flow equations to be derived later, the 
figure shows how these dimensionless couplings change 
under a change of the scale k. 

The interpretation of such a RG-flow diagram is merely 
textbook knowledge. We sketch an illustration in Fig. [2] 
to extract the essential features from Fig. [T] If the flow 
goes toward the left side to smaller values of ft 2 ,, the 
symmetric state at (p = becomes more and more unsta- 
ble and the system falls into the symmetry-broken phase 
with ip 0. If the flow goes in the opposite direction to- 
ward the right to larger values of on the other hand, 
the symmetric phase is a stable ground state. Therefore 
there is a critical line separating these two regimes. If the 
flow starts on the critical line, the system approaches the 
IR fixed-point which is the critical point characterizing 
the phase transition. 

In the shaded region in Fig. [5] the flow goes deeper 
into the region with p,Q > and Aio < 0. Usually such 
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FIG. 1. RG flow diagram in the fij,-\ik plane for the d = 
3 0(8) scalar theory that is realized from the U(2) x U(2) 
theory with A2A = chosen initially. The arrows indicate the 
direction from larger to smaller k and the dots represents the 
fixed points. 




FIG. 2. General structure of the RG flow for the d = 3 O(N) 
scalar theory. In the shaded region the flow goes to > 
and Aio < that is the "necessary condition" for a double- 
well shape in the effective potential JT]). 



flow patterns are physically meaningless in the analysis 
of the O(N) theory because negative A 10 makes the po- 
tential (TT]) unbounded and the theory is not well defined 
there. Nevertheless, from the interest in the first-order 
transition, this shaded region is interesting. This is be- 
cause, in view of the potential form ([1]), the destination 
of the flow, p,Q > and Aio < 0, is just the "necessary 
condition" to realize a double- well shape (see Fig. [3]). Of 
course this is not a "sufficient condition" and the ex- 
istence of the double-well shape depends on how large 
the sixth-order coupling constant £i associated with tp 3 
is (see Eq. (|35]l). If fi is too small, the potential is not 
stabilized until very large values of tp have been reached, 
for which an expansion like Eq. ([T]) is unsuitable. If £i is 
too large, the second minimum at tp 7^ simply vanishes 
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FIG. 3. Potential shape for fi k > and Aifc < 0. Depending 
on the sixth-order coefficient the potential may take a double- 
well form. 




and the symmetric phase is always the most stable. Nev- 
ertheless, we can say that, if there is a region where the 
first-order phase transition exists, the flow must satisfy 
the necessary condition, meaning that the flow must be 
headed for the region with Jxq > and Aio < 0. 

The U(2) x U(2) theory has at least one more coupling 
constant A2 in addition to /j, 2 and Ai (and (1) as defined 
in Sec. IIV1 For A2 = the theory is reduced into the 
0(8) scalar theory. Then, recalling the discussions in the 
previous paragraph, we can formulate our expectations 
about the evolution of the phase diagram. 

An interesting question is the following: How does the 
RG-flow diagram or more specifically the shaded region 
in Fig. [2] change in the U(2) x U(2) model as the addi- 
tional coupling A2 increases from zero? This question 
does not seem to have been answered explicitly in the 
earlier works by Berges, Tctradis, and Wetterich [42|, HU 
nor by Alford and March- Russell , though those works 
contain extensive and detailed discussions. 

Answering this question is quite intriguing because, as 
we have explained above, the shaded region has a direct 
connection to the emergence of the double-well shape in 
the effective potential. We can understand the results in 
a way parallel to the 0(8) case, using the phase diagram 
in space spanned by ft 2 and Ai^. Primarily this is just 
a replacement of the variables used to draw the phase 
diagram, but it also provides a useful interpretation in 
terms of the shape of the potential and we believe that 
it is worth taking a closer look at this interpretation in 
the following. 

Before we discuss the results shown in Figs. [4] [5j and 
[SI we remark that the horizontal and vertical axes are 
not, strictly speaking, identical to those in the previous 
Fig. [TJ In the diagram we use /2 A and Aia, which are 
the coupling parameters at the UV scale k = A, to show 
the shaded regions. If the flow starts with initial /2 A 
and Aia inside the shaded region and the specified A2A 
for each figure, it goes toward /2g > and Ai < at 
the IR scale k « 0. In the numerical calculation, we do 
not actually take the limit of k — > but stop the flow 
at fc/A = 0.1. This is a prescription for the numerical 
calculation (i.e. the choice of 0.1 is just arbitrary). As 



FIG. 4. Initial condition region at the UV scale necessary 
for the double-well shape in the effective potential at the IR 
scale. Results are from the fourth-order Taylor expansion. 
The line labelled with "lst-order" represents the phase tran- 
sition points of fluctuation-induced first order as a result of 
full numerical evaluation in the grid method. 



discussed in 0, [42], |43| . in the context of a first-order 
transition with bubble formation, it would make physical 
sense to stop the flow at some point whose scale should 
be related to the bubble thickness. 

The RG-flow diagram for the U(2) x U(2) model sim- 
ilar to Fig. [1] is available from the solution of the RG 
equations for the coupling parameters, that is, the Tay- 
lor expansion coefficients. In the same way as for the 
O(N) model we have performed the Taylor expansion for 
the U(2) x U(2) system to derive a set of partial differ- 
ential equations for ft 2 ,, Aifc, and A2A;. Then, solving the 
differential equations numerically, we have identified the 
shaded region corresponding to that in Fig.[5]in the O(A^) 
theory. 

Figure U shows the results from the Taylor expansion 
up to the fourth order in terms of the field (or the second 
order in terms of ip — <fi 2 ). It is clear that our expectation 
is manifestly fulfilled: The shaded region is elongated 
from the unphysical region in the O(N) case and spreads 
almost along the critical line toward the physical region 
at Ji 2 A < and Aia > 0. 

As we will describe in detail in a later section, we have 
carried out a full numerical calculation without resorting 
to an expansion by putting the function Uk(y>) on a grid 
of discretized field values. In this way we can confirm 
the double-well potential form already reported in the 
earlier works by Berges, Tetradis, and Wetterich [42ll43j. 
While our method does not surpass the earlier works in 
precision, we use a regulator that leads to much simpler 
flow equations. 

The line labelled with "lst-order" in Fig. @] represents 
the phase transition points of fluctuation-induced first 
order. The strength of the first-order transition weakens 
with decreasing /2 A and eventually it becomes smaller 
than the resolution in the grid method. We indicate this 
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FIG. 5. The shaded region indicates the initial conditions 
at the UV scale necessary to obtain the double-well shape in 
the effective potential at the IR scale. Results are from the 
sixth-order Taylor expansion for A2A = 0.4. 



FIG. 6. The shaded region indicates the initial conditions 
at the UV scale necessary to obtain the double-well shape in 
the effective potential at the IR scale. Results are from the 
sixth-order Taylor expansion for A2A = 0.7. 



III. FORMALISM 



by changing the solid line to a dotted one. We can see 
that the agreement between the results from the fourth- 
order Taylor expansion and from the full numerical eval- 
uation is limited to a qualitative level. Guided by the 
expectation that the agreement would be better, we have 
investigated the phase diagram from the sixth-order Tay- 
lor expansion. 

We show the results from the sixth-order Taylor ex- 
pansion in Fig. [5j The phase transition line seems to 
be smoothly connected to the shaded region in this case. 
The agreement is, however, not as good as we expected. 
[We have checked that the second minimum in the po- 
tential lies in the region where if is sufficiently smaller 
than unity in all the cases, so that the Taylor expansion 
should work in principle.] We will discuss the quantita- 
tive comparison for the shape of the potential with and 
without expansion later. 

So far, we have seen the results for a fixed choice of 
A2A = 0.4. According to [H, EH, the first-order phase 
transition is strengthened with increasing A2A. We can 
reconfirm this observation by repeating the calculations 
for larger A2A- Figure H] shows the results for a larger 
choice of A2A = 0.7. We can see that the shaded re- 
gion in Fig. [5] becomes slightly wider than that in Fig. [5] 
Also the dotted line disappears since the first-order phase 
transition becomes more prominent and is resolved by the 
grid in the region shown. 

If we choose even larger values of A2A, we find that the 
shaded region is pulled not along the critical line, but 
upward toward larger values of p,\, while the position of 
the first-order transition line hardly moves. 



To make our discussion as self-contained as possible, 
we shall briefly review the functional RG formalism and 
the equations we are using in this paper. The functional 
RG is an exact formulation of quantum field theories in 
the form of a functional differential equation [44|-[46j. In 
the formulation we are using, the central object is the 
effective action which is calculated non-perturbatively. In 
the functional RG one defines an average effective action 
with the IR-cutoff scale k. For a general scalar field \ 
whose classical counterpart is denoted by <fi = (x) here, 
the fc-dependent effective action is defined by 



log 



/ 



[dx] exp(-S[x] 



8<f>(x) 



AS k [ X 
[x(x) - 4>{x)] 



(2) 



where the IR-cutoff dependence is introduced by a mass- 
like term, 



AS k 



d d q 
(2n) d 



R k {q)tt>{-q)tt>{q) 



(3) 



in momentum space. Here we consider d-dimensionaJ 
space-time generally. We note that the momentum- 
dependent mass R k {q) should cut off IR modes with 
q < k. That is, R k (q) must satisfy the following require- 
ment (44H46| ; soft modes with g<i should be quenched 
by large R k (q) ~ k 2 , while R k (q) ~ for q k. Equa- 
tion ([2) together with the properties of R k (q) leads to 
the boundary condition for the average effective action; 
Ta [<f>] = S [<fi] where A is the UV scale where the RG flow 
starts, and Tq[4>] = T[4>] where r[(j>] is the full effective 
action defined in the textbook manner. In principle, the 
IR-cutoff function can be arbitrary as long as it satisfies 
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the above-mentioned properties. In this paper we make 
use of the "optimized" choice of Rk (<?) [H, [zH, IS pro- 
posed by Litim [74| , 



Rk(q) = (k 2 ~q 2 )0{k 2 



(4) 



where q refers to a d-dimensional vector in Euclidean 
space-time and 8 is the Heaviside step function. The 
choice of a cutoff function is in general not unique, but 
the optimized choice [z3] ensures the correct universal 
behavior (7j|. (For finite-T studies it is a better choice 
to use the cutoff for only spatial momenta [76l l77j|.) 

Armed with these definitions, we can reach the func- 
tional RG equation (known as the Wetterich equa- 
tion [44|) by differentiating Tk[4>] with respect to k, i.e. 



i. 



dk 



-Tr 



dRk 
dk 



(5) 



Here Tr means the trace over any intrinsic indices of <fi as 
well as the trace over the functional space. In the above 

(21 

; [4>] denotes the second-order functional derivative of 

(2) 

Tk [4>] with respect to the field cj>, and T], ' [(f)] + Rk is thus 
an inverse of the full propagator including the IR-cutoff 
term. From this one can acquire an intuitive interpreta- 
tion for the left-hand side of Eq. ([5]) as a one-loop inte- 
gration attached to the cutoff derivative, dRk/dk. Nev- 
ertheless, since the derivation assumes no approximation, 
Eq. ([5]) is "exact" and contains full quantum effects. 

In practice, however, solving the functional RG equa- 
tion without any truncation is as difficult as finding an 
exact solution of the theory. We should thus make an 
approximation. Here we adopt the local potential ap- 
proximation (LPA) in which the effective action Tk [4>] is 
assumed to be of the form, 



[0]= / dfixf^drf + UkW) 



(6) 



in d-dimensional Euclidean space-time. The fe-dependent 
effective potential Uk(4>) is a local function of the fields. 
This approach is frequently taken in theories involving 
only scalar fields and is sufficient as long as the anoma- 
lous dimension remains small. Thanks to the particular 
choice (U]) one can easily perform the ^-integration after 
Eq. © is substituted for Tk[4>] in Eq. (|5|) (and <j> is as- 
sumed to be spatially constant). Finally the functional 
RG equation simply reads 



dUk 

dk 



K, 



,k d+1 J2 



1 



(7) 



The index i refers to the field components fc. Here, we 
have defined — Sd/(2Tr) d with the volume of a d- 
dimensional sphere Sd given by 



S d 



2n d / 2 
dT{d/2) 



(8) 



that is, for example, for d = 3 we need 

f 



6tt 2 



(9) 



The energy in the denominator of Eq. is 

E 2 = k 2 + M 2 , (10) 

where M, denotes the mass eigenvalues of the second- 
derivative matrices of the effective potential, 



Mi, 



d 2 U k &) 
dcj) t d<f>j 
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Although Eq. Q is very simple, it encompasses rich con- 
tents of the critical phenomena, as we will see later. 



IV. U(2) x U(2) SCALAR THEORY 

We are now ready to proceed to our main subject, i.e. 
the fluctuation-induced first-order phase transition. For 
this purpose we adapt the U(2) x U(2) scalar theory in 
d = 3 dimensions and introduce some notations accord- 
ing to hadron physics. 



A. Notations 

We can write the Lagrangian density of the 
U(2) x U(2) theory conveniently as Q 



C = ^tr[9 M $^$t] + ^ 2 tr[$$t] 
- 5l (tr[$$ t ]) 2 - 5 2 tr 



(12) 



with a complex 2x2 matrix $. This Lagrangian density 
is invariant under the transformation, 



$ — ► V L <PV^ , 



(13) 



where Vl and Vjj are independent U(2) matrices. In the 
context of chiral symmetry of QCD, the matrix $ has 
a parametrization in terms of the hadronic degrees of 
freedom: 



S + ill = ^2 t a (<7 a + i^a) , 



(14) 



a=0 



where the t a are the u(2) generators (and unity) normal- 
ized according to tr = S a b- We can write 



7I C 



V2 



n 



7T + 71 77 



-75* +71 11 



(15) 
(16) 
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using the conventional notation in hadron physics; aP = 
a 3 , a ± = (a 1 =p ia 2 )/\/2 = (ai =p «cr 2 )/V2, er = cr for 



the scalar sector and ir ' 



7r± = (tt 1 T iir 2 )/V2 



(ir i TW2)/\/2, ?7 = To for the pseudo-scalar sector. [This 
77 does not correspond to the physical 77 in nature but to 
the flavor-singlet 770 which mixes partially with 773 to be- 
come the physical 77'.] With this notation the Lagrangian 
density in Euclidean space-time is expressed as 

C = d^ad^a + d^TT-d^TT + d^d^rj + d^a-d^a + UA , (17) 

where we hnd the potential, 

1 



Ua = 2A4(c + 7? + V + a ) 



+ [91+ o52)(^ 2 +7? 2 +7? 2 + a 2 ) 2 



(18) 



from Eq. l[T2| . To simplify the notation in what follows, 
we introduce the new couplings and variables; 



+ 2g 2 (cr 2 + 7t 2 )(t7 2 + a 2 ) - (0-77 + 7? • a) 2 , 



Ai = 4! (gi + -g 2 ^ , A 2 = 2g 2 , 
ip = cr 2 + 7? 2 + if + a 2 , 
C = (a 2 + tt 2 )(v 2 + a 2 ) - ((JT7 - tt • a) 2 



(19) 
(20) 



Roughly speaking cr plays the role of the component of 
4> that acquires a finite expectation value in the O(N) 
scalar theory and ip defined above is the counterpart of 
(f> 2 . Using these variables we can rewrite Eq. (fTB)) into a 
concise representation as follows: 

^A(^0 = ^AV+^AiA(y3 2 + A 2A C- (21) 



B. Functional Renormalization Group Equation 

The functional RG equation in this system with 
U(2) x U(2) symmetry has the same form as the generic 
one ([7]) in the previous section. Thus, the functional 
RG equation in the LPA takes the following form (with 
d = 3): 



dk 



(22) 



where the energies in the denominator arc for i = <x, 
if, 77, and a, given with the mass eigenvalues obtained 
from the potential curvature. One might imagine that 
this straightforward procedure should work if U k is given 
as an explicit function in terms of cr, if, 77, and a. In 
fact, however, the U(2) x U(2) symmetry constrains the 
potential and U k is a function of only two independent 
variables instead of four; that is, U k is parametrized by 
ip and £ only as exemplified in Eq. (|21[) . 

Below we will elucidate how to read off the expressions 
for the mass eigenvalues from the derivatives of U k (ip, £), 



which requires a particular strategy for the calculation. 
Actually taking the second derivative on Uk{<p,£) is just 
a simple manipulation. The most non-trivial part lies in 
the following: The resulting masses are easily expressed 
in terms of er, 7? , 77, and a, which must be converted into 
ip and £ in the end. 

We can simplify the problem a bit by observing that 
one field (cr in our case) acquires a non- vanishing expecta- 
tion value and others have a vanishing expectation value 
and can be set to zero in the end. This means that we 
should take ip — > a 2 and £ — > after writing down the 
functional RG equations. It is therefore a good trunca- 
tion to keepthe ^-dependent term only up to the first 
few orders 14211 : 



u k (<p,z) = v k (<p)+w k (<p)t+o(e 



(23) 



In the present case we truncate the above expansion up 
to the linear order in £ (where £ is by definition of quartic 
order with respect to the meson field variables), so that 
our results cover at a minimum the same behavior that 
is captured in the perturbative e-expansion. In other 
words, we can say that we have chosen an initial condition 
where the coefficient of the £ 2 -order term is zero (as in 
Eq. pT]) 1 ). and a truncation where it is set to remain 
zero throughout in the RG flow. This is the same-order 
truncation as the one used in [42|, |H| . 

Let us demonstrate our strategy by taking an example 
of calculating M 2 . The second derivative immediately 
leads to 



Mi 



d 2 (V k + W k Q 
da 2 

4a 2 (V k " + W^) + 2(Vl + WU) 
+ 8aW k [cr(?? 2 + a 2 ) - 77(0-7; - 



2W k a 2 
■a)] . 



(24) 



It is highly non-trivial how to rewrite this expression for 
M 2 into a form in terms of ip and £. We will do this by 
thinking of a particular situation where only a and a 1 
take a finite value. This is an almost unique choice for 
simplification; we can take 7? and 77 as small as we like, 
but a is constrained by Eq. (f20f for a given set of <p and 
£. Then we have 

p> = a 2 + (a 1 ) 2 , £ = aV) 2 , (25) 

from which wc can solve for a and a 1 as 



J- \2 



± 



which can be expanded in terms of £ into a form; 



a = p 



(a 1 ) 2 = ~ 



(26) 



(27) 



up to the linear order of £ (higher-order terms are unnec- 
essary in the £ — > limit). These expanded forms look 
singular at ip — > 0, but one should keep in mind that the 
expansion assumes tp 3> £ in Eq. (|26p and one can always 
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choose £ in such a way that £ is small enough to satisfy 
this inequality. Therefore the singularity at ip — > is 
only spurious, as we will check explicitly later. 
Substituting these expressions into M 2 we find 



Mi 



i 



(28) 



In the same way we can read the mass expressions off for 
other meson fields: 



= M 2 = 2VI + 2WU , 
M\ = 



M\+2wJ 
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Af Q 2 ! = Ml, + 4F fe "| + 8W££ + 2iy fe 
M 2 2 =M 2 3 =M 2 1+ 2W^-£) . 



(29) 
(30) 

(31) 

(32) 



We remark here that, if U k (tp, £) contains a term of order 
£ 2 , the mass M 2 1 would have an extra term proportional 
to <p£, which is vanishing in the present case due to our 
choice of the initial condition and the truncation. 

It is important to note that the mass matrix does not 
only have diagonal components but also off-diagonal ones 
which lead to a mixing between different fields. These 
off-diagonal components are given by 



Mi 



M 2 _ wl = 2W k jl 



(33) 
(34) 



In order to make use of the simple expression for the 
functional RG equation in Eq. (|22|) . we have to identify 
the eigenmodes (a, a 1 ) for the a-a 1 part and (jy,^ 1 ) for 
the rj-ir 1 part. It is not difficult to find the eigenvalues of 
the 2x2 matrix spanned in a-a 1 space and the results 
are found to be 



M~ = W k " ( tp - 4 ) + Wgrt + 2Vl + 10WU 



Ml 



2W k 



Z , 8(14" + w k + w^) 2 £ 



'r 



w k 



>- - 



w k + w kV f e 



2K" - W k 



(35) 



(36) 



In the same way, regarding the rj-ir 1 mixing, we diago- 
nalize the mass matrix, yielding 



M 2 



2Vi + 2W^ + 2W k Vl, (37) 
M 2 , = 2Vi + 2WU -2W k ^. (38) 



The remaining task is to decompose the functional RG 
equation into the one contributing to V k (tp) and the other 



contributing to W k (<p) . As seen from the LHS of Eq. (|22[) , 
the part contributing to W k {ip) must come from terms 
proportional to £. (In the above the terms of 0(y/^) 
cancel between M| and M? ± .) Thus, we should expand 
the RHS of Eq. (|22|) and identify the flow equations for 
V k ((p) and W k (ip) from the contributions of O(£ ) and 
0(^ x ), respectively. 

The 0(£°) contribution is easy to derive from Eq. 
that is, 



dVk 
dk 

where we define 



E 2 

7T7] 



E 2 ) 



El v = k 

E 2 = e 



k 2 + 2Vi + WZ'tp 



+■ , 

2^ 



2W k <p 



(39) 



(40) 
(41) 
(42) 



The calculations in the order 0(^ 1 ) are rather compli- 
cated. After tedious but straightforward computations 
we finally obtain 



dW k 
dk 



-Kt k 



d+l 



8W> 



■4W k ip~ 



E 4 



SW 2 
E e 



-4V fc 'V^ 1 + 4WfeV + lOWfc + 2W k tp~ 1 



8(V fe " + iy fc + W» 2 x 



2^' 



W k 



-Lp 



4V fe 'V" 1 + 14Wfc - QWkip' 1 



W + w * + 
2V" - W k 



1 

Ei 



-El}- (43) 



The above equations ([5^1 . and (|4*3^) are our cen- 

tral results at the algebraic level. We emphasize that, 
thanks to the choice of the optimized IR regulator, these 
expressions no longer contain momentum integrals and 
are much simpler compared to those in earlier works. 
This also makes the connection to perturbative RG re- 
sults more apparent. 

Now let us make sure that Eq. (|4"3")l is not singular at 
ip — > 0. We see that the energy denominator becomes the 
same in all terms, i.e. E 2 = E 2 V = E 2 , as the limit p — > 
is taken. Then, we pick up the singular terms from the 
parentheses in the RHS of Eq. (|4"3")l . whose coefficients 
amount to 



iW k + 



-4V fe " + 2W k 



Kv^ + w k f 



6W k 



W k 

(v k + w k ) 



- 1 



W k " 



W k 



(44) 



0, 



where the first term is due to the part with £^L, the 
second parenthesis to the part with E 4 , and the last 
parenthesis to the part with E 4 in the functional RG 
equation (|4"3")l . 
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C. Flow in Parameter Space 



Let us check that our equations (|39|) , (|42[) . and 
(|4"5|) arc consistent with the known results from the e- 
expansion H, 0] ■ To this end we expand the potential up 
to the fourth order or sixth order in terms of the fields: 



X 2 k + &k<P , 

(45) 

where a constant offset of the potential energy is dropped. 
It should be noted that £ is already of fourth order. For 
notational convenience, we will use Aifc and \2k and will 
change to the couplings g\ and g 2 only at the end to make 
a comparison with the e-expansion results. 

By expanding the flow equation in terms of ip, we can 
formulate the flow of the coupling constants (where we 
will set Cifc = C,2k = and work up to the fourth order 
for the moment) as 

dVk 
dk 



(46) 




4X 2 k + 18X lk X 



2k 



9E e 



where E 2 — k 2 + [i 2 . Wc also have 

dW k 8X 2k d+1 8X lk X 2k + l6X 2 k 

-W = -Bk =Kih m • (47) 

We note that there is a cancellation which eliminates 
the terms proportional to (p^ 1 from the results, as wc 
have already checked before. This allows us to easily 
decompose the flow equations into 



dfxj _ 2K d k d+1 
dk ~ 2,E A 

dX lk = 16K d k d+1 
dk 3E 6 ' 

dX 2k 8K d k d+1 



(5A 



lk 



18A 



2 k 



(48) 



'2X 2 k + 9X lk X 2k + 54A 2 fe ) , (49) 



(AifeA 2 fe + 2X^ k ) 



(50) 



dk E 6 

up to the fourth order for d dimensions. At this point, it 
is a simple task to convert the above into the e-expansion 
results obtained by Pisarski and Wilczek Q. We rescale 
the coupling constants in the following way (see also 
Eq. (OH): 



Ai fc = 4!y(.9 1 + i5 2 )fc 4 - d , 



hk = ^92k 4 - d , (51) 

and use £l 2 = k~ 2 f-i 2 as previously defined. Then, after 
a short calculation, in d = 4 — s dimensions, we readily 
arrive at 



k dk 

fc ~Hr = -f ^ 2 + 2 ^ 2 + \® 2 ' 



8 8 

-m + 3<7i + 3.91.92 + al 



(52) 
(53) 



up to quadratic order in e on the RHS, where we eval- 
uated K d at the expansion point d = 4 and used K 4 = 
(327T 2 )™ 1 . This clearly shows that our calculations cor- 
rectly reproduce the results in the e-expansion in the 
leading order, as they should. One can prove that the 
flow with respect to g\ and §2 (with \x k = 0) does not 
have any stable IR fixed point |7|. This is a reasonable 
but indirect argument that the phase transition should 
be a fluctuation-induced first order one 0. It is highly 
non-trivial why the behavior of g~i and g 2 is capable of 
describing a form of the potential which admits a first- 
order phase transition. Let us consider this question by 
taking account of the flow behavior of p, 2 , which already 
goes beyond the e-expansion (in which /2| = 0). 

The set of equations (|48]), (J49J), and ([50]) is our con- 
venient starting point to deal with the parameter flow 
because one can intuitively connect Ai^ and A2fe to the 
shape of the potential; Aifc is a coefficient of the quadratic 
term ip 2 (i.e. the quartic term c 4 ). The flow equations 
for the dimensionless couplings in d = 3 dimensions, 
Aifc, and A 2 fc, are 



fc M - -OH 2 - 

k dk - 

dXi k 



1 



(5Ai fc + 18A 2fe ) , (54) 



9^(1+^)3 
k^r^- = ~X 2k + 



9^(1+^)2 

(2X 2 k + 9X lk X 2k + 54X 2 2k ) , (55) 



dk 



37r 2 (l + /2 2 fc ) 3 



;A lfc A 2fc + 2A 2 fc ) . (56) 



In the same way it is tedious but straightforward to write 
down the differential equations including sixth-order cou- 
plings Cifc and (2/i. The expressions are easily derived 
from Eqs. (|3T))) and (|4"3")l and they are too lengthy to dis- 
play here, so we shall give their explicit forms in the 
Appendix. 

We solve these differential equations numerically. The 
RHS of these equations has no explicit dependence on the 
scale k and the LHS can be rewritten as kd/dk = d/dt 
using t = ln(fc/A), so that we can describe the RG flow in 
terms of this variable. Wc numerically trace the running 
of these couplings from t — to t — —2.3 (i.e. k = 0.1A) 
to determine the flow diagram, using the fourth-order 
Runge-Kutta method. We used an adaptive step size 
(< 10~ 4 ) so that no numerical instability occurs. [In 
later discussions we will use the fifth-order Runge-Kutta- 
Fehlberg (4,5) (RKF45) method to solve the functional 
equations. In the present case, where the partial differ- 
ential equations are relatively simple, the fourth-order 
method is sufficient.] 

The case with X 2 \ = (Fig. [JJ); Let us begin with a 
simple case. We see that A2A; = is a solution of Eq. (|56"j) . 
If we start the RG flow with the initial condition A2A = 0, 
this means that X 2k remains zero for any k throughout 
the flow and Eqs. (|5~4"|) and (|55|) determine the evolution 
of ft 2 and Aifc with changing k (up to the fourth-order of 
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the Taylor expansion). In this case with \ 2 k = the flow 
equations are reduced to those in the 0(8)-symmetric 
scalar theory. It is easy to make sure that we can also re- 
produce the well-known results from the e-expansion by 
setting X2k = and using K4 = (327r 2 ) -1 in the above 
expressions. To go to the next-to-leading order in the 
e-expansion one must go beyond the local potential ap- 
proximation, as done in the derivative expansion [461 UM 
or in the so-called BMW approximation (79|, which is 
beyond the scope of our present approach. 

Using Eqs. (JSJ) and ([52]) with \ 2k = 0, we plot the RG 
flow to obtain Fig. [TJ We find an IR fixed point at 



5 

37 



-0.135 



— 9tt 2 /32\3 



(57) 



We find indeed in Fig. Q] a fixed point corresponding to 
the second-order phase transition at these values of the 
couplings. It is obvious from Fig. [T] that there is no RG 
trajectory flowing from the (/2 2 < 0, Aifc > 0) region 
into the (/z 2 . > 0, Aifc < 0) region. Thus it is not possible 
to arrive at a fluctuation-induced first-order phase transi- 
tion as long as the initial A2A is chosen to be zero. We see 
that such a flow which would favor a fluctuation-induced 
first-order transition is prevented by the presence of the 
UV fixed point at the origin from which the flow emerges. 



The case with A2A 7^ (Figs. [5l Now we pro- 
ceed to the case with a finite A2fc. Interestingly enough, 
the flow pattern changes drastically once A2A deviates 
from zero. Then we have to analyze the flow pattern as 
a function of /I 2 ,, Aifc, and X 2 k- To visualize the flow pat- 
tern, we solve Eqs. (|54[) . (|55p . and (|56[) for a particular 
initial condition with a certain choice of A2A, and inves- 
tigate the destination of the flow in the /1 2 -Aifc plane. As 
we have already elucidated in Fig. [5]in the 0(8) case and 
also in Figs. 01 [5J and [6l we can classify the parameter 
regions according to qualitative features of the potential 
form; 



p? k < and Au > 
p, 2 k > and Aifc > 
fil > and Aifc < 



Symmetry-broken 
Symmetric 

Necessary for Double-well 



as k — > 0. 

We have already discussed our results in Sec. HT1 As a 
final remark in this subsection, we point out an interest- 
ing consequence from the flow equations. 

From an analysis of the flow equations (|54"|) , ((53]) , and 
(|56p we can indeed confirm that in this approximation no 
fixed point with A2fc 7^ exists in d = 3 dimensions, and 
in the presence of a non-zero coupling A2fc a second-order 
phase transition associated with this fixed point cannot 
be accommodated. 

To see this, the fixed point (/2 2 , Ai*, A 2 *) is located by 
imposing the condition that the RHS's of Eqs. (|54]l. ([55]) . 
and (|56p vanish. For A2* 7^ 0, the fixed point condition 



on 



-0.0001 




FIG. 7. Evolution of the form of the potential with decreasing 
k/A in the case with A2A = 0.4. The solid curves represent 
the results from our grid method (see the text for details) 
and the squares represent the results obtained in the method 
formulated in [sjj. The potential is given as a function of 
^/7p. The initial parameters are chosen as /i A = —0.05 and 
Aia = 0.35. 



can be simplified and the resulting system of three equa- 
tions for the three couplings /I 2 , Ai*, and A2* can be 
solved analytically. We find that for d = 3 no solutions 
for real values of the couplings exist (that is, /x 2 , Ai*, and 
A2* are all complex then) , and that hence the RG flow in 
this approximation docs not admit a second-order phase 
transition. 



D. Functional Solution 

Without an evaluation of the higher-order terms we 
cannot locate exactly where the fluctuation-induced first- 
order phase transition takes place on Figs. 01 [5l and |6] 
Nevertheless, it is clear that the second-order critical line 
is overridden by the region which satisfies the condition 
for a double-well potential as the initial value for A2A 
grows (see particularly Figs. Q] and 0}. 

We can check if the analytical study in the previous 
subsection is qualitatively correct by looking at the full 
functional solution obtained with the grid method. Since 
there is no scaling property expected in the case of the 
fluctuation-induced first-order phase transition, we plot 
the results in terms of unsealed variables (made dimen- 
sionless not with k but by the UV scale A). We give the 
numerical values for /i 2 , Aifc, and A2fc in units of A and 
we omit A hereafter. 

In the numerical calculations, we used a simple for- 
mulation for the grid method. We first discretize the 
function Uk(y>) in terms of a = ^JTp instead of ip, since 
we are interested in the characteristic form of the effec- 
tive potential in the vicinity of tp ~ or a ~ 0. It is 
thus more suitable to use a to attain more resolution in 
the small-y region. Our mesh size is Act = 10~ 2 and the 
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mesh spans the interval between <j = ±2. The singular- 
ity in the flow equations at a = is only spurious due to 
cancellation as we have checked. In the numerical calcu- 
lation one should remove this singularity analytically or 
discretize the fields not to hit the origin. We here took 
the latter, i.e. we choose the support points such that 
the origin is not included; <7j = Aa(i + 0.5) with i rang- 
ing from —200 to 199. To evaluate the derivatives U' k 
and U k we used the 7-point formula. At the boundaries, 
where we cannot take 7 points (i.e. 3 adjacent sites from 
the point where the derivatives are evaluated), we used 
the 5-point formula and the 3-point formula, and eventu- 
ally at the very end we approximately used the next-site 
value. These procedures enhance numerical errors locally 
near the edges, but if the edges are sufficiently far from 
the small- (p region of our interest, those errors do not 
influence the results. The step size Ak is an adaptive 
variable so that the computation proceeds without nu- 
merical instability. Typically Afc is less than 10~ 4 . Then 
we found no instability in the fourth-order Runge-Kutta 
method. In summary, after discretization we treat the 
function Uk(f) as an array labelled by site number i and 
then proceed as follows: 

1. Choose an initial value of k = A and fix the po- 
tential parameters /ia, Aia, and A2A which also de- 
termines the initial configuration of the array for 
Va(ct) and Wa(ct) (where a = ^fip). 

2. Evaluate the RHS of Eqs. (J39J) and (|43j) using V k (<r) 
and Wk(cr), and the numerical derivatives. Up- 
date the potential from Vk(cr) and Wk(cr) down 
to Vk-Ak(cr) and Wk~Ak{cr) with the Runge-Kutta 
method. 

3. Iterate the above calculation until the scale k 
reaches zero or a sufficiently small value, so that 
the full effective potential Vb(cr) and Wo(cr) results 
at the end from the functional RG flow. 

Such a discretization scheme may look very simple 
compared to more sophisticated processes such as the 
one proposed in 18011 , which has actually been used in 
the earlier works [42|43|]. In this method Uj~ and U' k are 
solved for by the evolution equations and U k and U' k " are 
determined by the matching condition. While the idea 
sounds very different from the simple numerical deriva- 
tives that we used, the matching condition leads in effect 
to a natural generalization of the formula for the numeri- 
cal differentiation. The accuracy is likely better than that 
of the the 7-point formula. In addition, for the method 
of H3|, usually the RKF45 method (fifth-order Runge- 
Kutta method with the fourth-order embedded for the 
error control) is used. This sophisticated method is very 
powerful, but at the same time, it is difficult to apply it to 
solve a flow equation with many terms, such as Eq. (|43| . 

Wc prefer to use the former simpler method as long 
as it works and we do not need high precision data. In 
fact we have solved the RG equations making use of both 
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V ip 

FIG. 8. Evolution of the form of the potential in the region 
with k/A < 0.2 in the case with A2A = 0.4 as k/A decreases. 
At k/A = 1 the initial parameters are chosen as fi\ = —0.05 
and Aia = 0.35. 

the simple and sophisticated method, and we present a 
quantitative comparison in Fig. [7] We clearly see that the 
potential exhibits a double-well form as a result of fluc- 
tuations and the deviation by two methods is not visible 
in the figure. We can conclude that our method is good 
enough for our purpose. 

We have already discussed our central results in Sec.HII 
and so we do not reiterate the discussion here. For the 
rest of this subsection we mention on some notable fea- 
tures in the potential evolution. 

Figures [7] and [5] show the evolution of the shape of the 
potential in the case of A2A = 0.4 starting with the initial 
condition /i A = —0.05 and Aia = 0.35, which is almost 
at the first-order phase transition point (as indicated by 
the solid line in Figs. 0] and [5]). We stopped the evolution 
at k/A = 0.1, which is legitimate because the effective 
potential becomes convex when k/A goes to zero and the 
double- well shape of the potential is less manifest [8lj . 
Here this happens below k/A ~ 0.1, as shown in Fig. [5] 
We can clearly see in Fig. [8] that the location of the min- 
ima does not change for k/A < 0.1, although the bump 
of the potential is becoming flat with smaller k/A. This 
is completely consistent with the observation in Ref. [42| 
and justifies our prescription to stop the RG evolution at 
k/ A = 0.1 in order to judge if the potential has a double- 
well shape. Even though we choose the step size Afc as 
an adaptive variable, moreover, our algorithm loses sta- 
bility at k/A ~ 0.02 and cannot reach the completely 
flat shape. We observed the same behavior even using 
the matching method with the RKF45 method. Prob- 
ably the implicit method (i.e. Lax method) might cure 
this problem. In any case, since our aim is simply to de- 
termine the transition line on the phase diagram shown 
in Figs. [?]and [5l there is no practical need to require the 
potential form at k/A < 0.1. 

The first-order phase transition occurs with varying 
/i A and/or Aia for a given A2A- Let us fix ji\ at —0.05 
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FIG. 9. Change of the form of the potential as the initial 
value Aia changes around 0.35. All results for the potential 
are given at fixed k/A = 0.1. Other initial parameters are all 
chosen to be the same as in Fig. 



present Fig. [10] to show the expanded results according 
to Eq. (|4"5|) together with the outputs from the grid cal- 
culation. It may be surprising at first glance that the 
curvature at ^JTp ~ has such large deviations from the 
fourth-order results and the sixth-order or full grid re- 
sults. We observe that the agreement in the curvature 
stays good as long as k/A is not too small (k/A ~ 0.5 for 
instance), for which the curvature is still negative. Even 
though the Taylor expansion in terms of the field ip should 
in principle work well in the small-<^> region, the fourth- 
order expansion is not good enough to capture the first- 
order phase transition quantitatively. We can see that 
the agreement in the curvature is significantly improved 
by proceeding to the sixth-order expansion, though the 
agreement in the quartic coefficient (i.e. Aifc) seems to be 
poor, which might be improved with inclusion of further 
higher-order contributions. 



V. CONCLUSIONS 



0.0001 
0.00005 


S? -0.00005 





6th-order ..♦** 


k/A =0.2 






4th-order 




k/A=l \ 











0.05 0.1 0.15 0.2 0.25 0.3 



FIG. 10. Comparison between the full effective potentials at 
k/A = 0.2 and the expanded results from the Taylor expan- 
sion up to fourth order and sixth order. The effective potential 
at k/A = 1 is presented for reference. 

here and change Aia to see how the first-order transition 
arises. Our results are depicted in Fig. [9] Again, we have 
stopped the evolution at k/A = 0.1, which is sufficient 
to find a first-order phase transition. It is obvious that 
there is a peculiar change of the shape of the potential 
which is associated with the first-order phase transition 
from Aia = 0.34 to 0.36. More precisely, we can locate 
the transition point by detecting the second minimum 
height and using the midpoint method for Aia, which 
leads to a critical value of Aia = 0.3509. We picked just 
one example of particular values of /i\ but can repeat 
the same procedure with different parameters to find a 
first-order phase transition line as indicated in Figs.|U[5j 
and [5] 

Finally let us discuss the comparison between the full 
functional solutions and the Taylor expansion results. We 



We have applied the functional RG technique in the 
form of Wetterich's flow equation to an analysis of a 
scalar theory with U(2) x U(2) symmetry. We have been 
able to recover the perturbative results by means of ex- 
pansion in terms of the field, which also encompasses the 
standard e-expansion results Q ■ The choice of a suitable 
regulator function and the resulting simple form of the 
RG flow equations make this connection apparent. Going 
beyond a perturbative expansion in small couplings, we 
have further used the functional RG flow equation in the 
local potential approximation to study the emergence of 
a first-order phase transition. To this end, we have dis- 
cretized the effective potential on a mesh grid to obtain 
the global shape of the potential for a set of choices for 
the initial conditions. Our results confirm those of earlier 
studies [Hm. 

In the U(2) x U(2) scalar theory there are two fourth- 
order couplings - Aife appearing in the quartic coefficient 
in the effective potential and A2fe which has no counter- 
part in the O(N) theory - in addition to the curvature 
fif,. We have found that the second-order critical line 
in the fil-Xik plane is gradually overridden by the first- 
order phase transition as the initial value A2A starts to 
differ from zero. We find analytically that no fixed point 
associated with a second-order transition exists for d = 3 
and A2fc 7^ in our approximation which is perturbative 
but beyond the e-expansion. 

In view of the expanded potential, a double- well shape 
of the potential is possible in the region /2 A > and 
Aia < 0, which would admit a first-order transition if 
the potential is stabilized by higher-order terms, which 
we investigated up to the sixth-order expansion. We note 
that a similar analysis has been performed in Ref. (Io| up 
to the eighth-order expansion for the Coleman- Weinberg 
potential and our results are qualitatively consistent with 
those of Ref. [l(| • From the flow diagram we learn that 
one way to understand the fluctuation-induced first-order 
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phase transition is that this first-order phase transition 
at the tree level with Ji\ > and Aia < penetrates into 
other (physical) parameter regions for A 2 a ^ 0. This is 
an interesting point of view and gives an intuitive insight 
into the emergence of a fluctuation-induced first-order 
phase transition. 

Regarding the application to QCD physics, it would be 
intriguing to study the chiral phase transition with both 
the meson fluctuations [H, [70] and the effective restora- 
tion of U(1)a symmetry taken into account. So far, for 
example, the so-called Columbia diagram in the mass 
plane of light and heavy quarks has been investigated as 
a function of the U(l)A-breaking strength [64], |62|. The 
possibility of the fluctuation-induced first-order phase 
transition, however, has not yet been considered because 
the model study was at the mean-field level. The present 
work gives a theoretical framework necessary for such in- 
vestigations in the future. 
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Appendix A: FLOW EQUATIONS IN THE 
SIXTH-ORDER TAYLOR EXPANSION 



Here we list a set of differential equations involving the 
sixth-order Taylor coefficients £i and £2- We note that 
Eq. (|4*5|) for /i 2 has no modification. Equations (|4"9"j) and 
(|5D1) are replaced, respectively, by 



dAi fc 16K d fc d +i 2 



2k 



dk 



8K d k d ^ 



54A 2 



2/v 



9AifcA2fc 
-27£ 2 (12Ci fc + C 2fe ) 
(A lfc A 2fe + 2A 2 fe - AC 2k E 2 ) , 



which arc reduced to Eqs. (|49[) and ([50)) for £ifc = (2k 
The differential equations for £1/. and (,2k are 



(Al) 
(A2) 

= 0. 





K d k d ^ 


dk 


108£ 8 


d( 2k 


4K d k d+ 


dk 


3E 8 



{I7\l k + 54A 2 fc A 2fc + 648A lfe A 2 fc + 2592A^ - 216£ 2 [2& fc (llAi fc + 18A 2fe ) + ( 2k (\ lk + 12A 2fc )] } , 

(A3) 



■{3A 2fe (2A 2 , + llA lfe A 2fe + 6A 2 fe ) - £ 2 [432Ci fc A 2fc + C 2fc (23A lfc + 114A 2fe )] } 



(A4) 



which were used to obtain the sixth-order results shown 
in Figs. [6] and OB 
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